Quench dynamics in spin crossover induced by high pressure 
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In this paper we have studied analytically and numerically dynamics of spin crossover induced 
by time-dependent pressure. We show that quasi static pressure, with a slow dependence on time, 
yields a spin crossover leading to transition from the state of quantum system with high spin (HS) 
to the low spin (LS). However, a quench dynamics under shock- wave load is more complicated. The 
final state of the system depends on the amplitude and pulse velocity, resulting in the mixture of 
the HS and LS states. 
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I. INTRODUCTION 

Spin crossover in condensed matter physics is a trans- 
formation of a system with one spin Si at each lattice 
site into another state with spin 62 induced by some ex- 
ternal parameter like strong magnetic field, high pres- 
sure etc. It accompanies by the energy level E\ and E% 
crossing, where E a is the local energy of the magnetic 
ion with spin S a (a = 1,2). Recently spin crossovers 
in magnetic oxides have been found under high pressure 
in FeBO^, CdFe 3 (B0 3 )^, BiFe0 3 3 , Fe 3 4 -- Below 
the Curie temperature of magnetic order spin crossover is 
accompanied by the sharp change of the magnetization, 
nevertheless it may be observed in the paramagnetic state 
like in CdFe 3 {B0 3 )^ as the sharp change of the XES 
satellite/main peak intensity ratio with pressure increase. 

Energy levels crossing results in the lost of analyticity 
in the energy spectrum at the critical point (in the ther- 
modynamic limit) 5 . Near the critical point adiabatic- 
ity breaks down and non-equilibrium phenomena asso- 
ciated with the drastically grown quantum fluctuations 
can drive the system away from the ground state. The 
final result depends on how fast the transition occurs. If 
the quench process is sufficiently fast, large numbers of 
topological defects are created and the final state, being 
characterized by mixture of high and low spin phases, can 
be essentially different from that been obtained as result 
of slow evolution. Qualitatively, quench dynamics can be 
described by the Kibble-Zurek theory of nonequilibrium 
phase transitions^—. 

In this paper we consider quench dynamics in spin 
crossover induced by time-dependent pressure. The pa- 
per is organized as follows. In Sec. II, a general model of 
spin crossover under high pressure is introduced. In Sec. 
Ill, we study quench dynamics. We consider two cases: 
a) Pressure is a linear function of time; b) Pressure is 
defined by a pulse of a given shape. We conclude in Sec. 



IV with a discussion of our results. 



II. MODEL 

The multielectron ion in a crystal field has the energies 
of terms for d n configurations determined numerically by 
the Tanabe-Sugano diagrams - as a solution of the eigen- 
value problem. Simple analytical calculations of the low 
energy terms with different spin value that is sufficient 
to study spin crossover has been done recentljJ^. The 
crystal field parameter increase linearly with pressure P. 
Thus the multielectron energies for spin S\ and S2 ( E\ 
and E2 ) are also linear functions of P. To distinguish two 
different spin states in the lattice we introduce the Ising 
pseudospin states \i) and | — i) for \df,S\) and \d™,S\), 
where i runs over all sites in the lattice. Thus we neglect 
the spin degeneracy of the d™ terms but capture the pos- 
sibility of energy level crossing that is the essential part 
of the spin crossover. Then, in the basis | + i), \ — i), the 
Hamiltonian of the system can be written as follows 



(1) 



where Aq = {E\+E l 2 )/2, e, = (E\-E\)/2, and 1, a z are 
the identity and Pauli matrices, respectively; the Hamil- 
tonian of interaction between the spins being Hij . 

The includes the isotropic Heisenberg term with 
the exchange interaction between nearest spins and 
the anisotropic term Ha- The interatomic interaction 
Iij is negligibly small in comparison with the interatomic 
Hund's coupling (ratio 10~ 2 ). Thus its contribution to 
the localized spin energy E\ and E2 due to the effective 
molecular field can be neglected. Nevertheless the ex- 
change interaction plays very important role: it results 
in the long range order and synchronize each spin in the 
same quantum state providing a cooperative behavior of 
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the spin system. If it were the ferromagnetic interaction, 
each spin at T=0 would have the maximal projection S 
with integer magnetic moment 2S. 

In all examples given above there is the antiferro- 
magnetic interaction. The ground state of the isotropic 
Heisenberg antiferromagnet has non integer local mag- 
netic moment due to the quantum spin fluctuations. It 
is known that for large spin S effect of quantum fluctua- 
tions is less important then for small spin, and for FeBOs 
spin is 5/2. Moreover the magnetic anisotropy addition- 
ally suppresses quantum fluctuations. For example, in 
FeBO-i the anisotropy field is 0. 3T— 1 and the measured 
value of the effective moment 2yJ S(S + I) = 5.9 is very 
close to the calculated for S = 5/2 value 5.916. 

Thus we conclude that due to anisotropy the magnetic 
moment at T = has integer value ( of course it is a 
property of the magnetic insulator that does not hold for 
itinerant magnets), and due to exchange interaction all 
spins are in the same quantum state. So spin crossover 
at T = is the transition of the whole crystal from one 
magnetically ordered state to another. Nevertheless the 
criterium of the transition can be found from consider- 
ation of the single ion energies crossover due to space 
uniform cooperative magnetic order. Anisotropic rela- 
tivistic interactions, for example a spin-orbital interac- 
tion, are also important because can mix different spin 
states inside single ion. 

We consider spin crossover far from the thermody- 
namic phase transition in the paramagnetic phase, it al- 
lows us to simplify this interaction and substitute the ef- 
fect of exchange with the effective mean field. This mean 
field is spatially uniform for the ferromagnetic insulator 
or two-sublattice for the antiferromagnet one. Examples 
given abovei~— correspond to the anti- or ferrimagnetics. 
In any case this mean field just renormalizes the interi- 
onic multielectron energies E\ and E%, and is irrelevant 
to the crossover phenomenon. Another interaction that 
is smaller then the exchange one is given by relativistic 
anisotropy contribution to the Hij. For example a spin- 
orbital interaction can mix different spin states inside 
single ion, and it occurs to be important in our problem. 

In what follows we will consider the simplified spa- 
tially uniform models. Motivation for this simplifica- 
tion is as follows. Despite that spin crossover is related 
to many body system, the essential features of its dy- 
namics can be described by effective, Landau-Zener type, 
Hamiltonian^. 

The Hamiltonian of our model is given by % = 
SiLi^i) where 

*-(o£) + U'C)- < 2 > 

The energy spectrum is given by e± = Ao ± a/e 2 + p 2 . 
Both Ao and e are pressure dependent. Further we as- 
sume that the spin excitation gap is given by 

E (P)= £o (l-0 (3) 



The crossover takes place at the point P c when e(P c ) = 0. 
The spin-orbit coupling A = pe lv (with p <C £o) mixes 
the different spin states, and it plays the role of quan- 
tum fluctuations in our Ising pseudospin basis. The di- 
mensionless spin gap, P/P c — f, plays a role of relative 
temperature (T/T c — 1) near the critical point T^. 

III. QUENCH DYNAMICS INDUCED BY 
HIGH PRESSURE 

We consider time dependent Schrodinger equation for 
the Hamiltonian ([2]) assuming for simplicity that the 
pressure is the linear function of time, P = P c (l + t/rq). 
Inserting this expression into Eq. @, we obtain e{t) = 
—Eot/rQ. The parameter tq depends on P and can be 
written as tq = P c /P. 

Let 1 1) and |0) are eigenstates of the operator <j z , so 
that CTz 1 1) = |1) and a z \0) — — 10). Expressing a generic 
state vector as 

\^{t)) = e-^ Ao(t)dt (Ci(t)e-^ /2 |I) + Co{t)e lv/2 \0)), 

(4) 

we find that the coefficients C\(t) and Co(t) satisfy the 
Schrodinger equation with the time-dependent Hamilto- 
nian in the Landau-Zener (LZ) form (in units h — 1) 

,d (C 1 {t)\_( -At P \ (Ci(t)\ . . 
l dt \C (t) J~{ P At) [C (t) w 

where A = Eq/tq. 

In terms of dimensionless scaled time r = y/At = 

t (P/P c - I) with t = y/WrQ = \J EoPc/P, the Landau- 
Zener model is described by the Hamiltonian 

«=(■;")■ < 6 > 

where ui = p/VA = tqp/eq is the dimensionless coupling 
constant. Writing \u(t)) = Ci(r)|l) +Co(r)|0), one can 
recast the Schrodinger equation ([5]) as 

i±\u( T ))=U(r)\u(r)). (7) 

Here the time r runs from the initial time Ti — —^EqTq, 
corresponding to the initial pressure Pi = 0, to final r/ = 
(Pf/P c — l)y/e^TQ, corresponding to Pf at the end of 
quench (Pf > P c ). Further we assume that £$tq S> I , 
then time r can be extended to ±oo, and the problem 
becomes fully equivalent to the LZ problem. 

The energy spectrum of the Hamiltonian (jSJ) is given by 
£±{t) = ±\/t 2 + ui 2 , and its instantaneous eigenvectors 
can be written as 
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FIG. 1: (Color online) Bloch vector's evolution (ui = 3). The 
motion starts at the south pole of the sphere and asymptoti- 
cally ends at the north pole. 



FIG. 2: (Color online) Bloch vector's dynamics (uj = 0.75). 
The evolution starts at the south pole of the sphere. 



where cos 0(t) 



UJ' 



The energy gap between 



the ground and excited states equals 2Vt 2 + uj 2 . 

From Eq. ([8]) it follows that while the ground state be- 
haves at r = ±oo as \u-(— oo)) — > |0) and |tt_(+oo)) — > 
|1), the excited state behaves as follows: oo)) — > |1) 

and |u+(+oo)) — > |0). The state oo)) corresponds 

to the high spin (HS) of the system and \u + (— oo)) cor- 
responds to the low spin (LS). Thus, if the system ini- 
tially was in the HS state, at the end of the evolution its 
ground state corresponds to the LS. Tunnelling between 
the positive and negative energy eigenstates, leading to 
the mixture of HS and LS, happens in the neighbourhood 
of the critical point r c = (P = P c ) when r g (— ui, ui)^. 

We assume further that the evolution of the system 
starts at the moment of time r, = —^/£o T Q (P( T i) = 0) 
from the ground state |u_(Tj)). Since p <C eq, we have 
|u_(Tj)) oc |0). This yields the following initial condi- 
tions: Co(rj) = 1 and Ci(tj) = 0. At the end of evolution 
we obtain, |n_(r) — > |1), while r — > +oo. 

In Figs. [T]-[3]we present the results of numerical solu- 
tion of the Schrodinger equation. In Fig. [1] [2] time evo- 
lution of the Bloch vector, n = (u\er\u), is shown. The 
motion begins at the south pole of the two-dimensional 
sphere S* 2 and for uj — 3 ends at the north pole. However 
for the choice of the parameter uj < 1 , numerical simula- 
tion shows that the Bloch vector never reaches the north 
pole, which corresponds to the LS state. This implies 
that at the end of evolution the quantum system does 
not remains in the ground state and its final state is the 
mixture of the LH and LS states. 

In Fig. [3] the probability P T = \Ci(t)\ 2 of transition 
|0) — > |1) is depicted for various values of uj. As can 
be seen, with decreasing of adiabadicity parameter uj the 
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FIG. 3: (Color online) Probability of transition P T obtained 
from exact solution of LZ problem as function of r. From up 
to down: ui = 3, 0.5, 0.25. 



transition probability decreases as well. Its asymptotic 
behaviour is described by the LZ formula (lli 



A. Exact solution of the Landau-Zener problem 

The exact solution of the Eq.© is given in terms of the 
parabolic cylinder functions^—, D_ i-i W 2/a(z), where 
z = v / 2 / re~ lTr / 4 . Assuming that initially the system was 
in the ground state, oo)) — > |0), we obtain the ini- 

tial conditions as follows: Cq(— oo) = 1 and C\(— oo) = 0. 
The probability of transition to the state 1 1) at time r is 
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given by-i£ — 

Pr = ye-™ 2 /4|^_ i _^ 2/2(r y2e 3 "/ 4 )| 2 . (9) 

As can be shown, the condition Eo ^ p may be recast 
to | Ti| 3> oj, where r, = — ^/£qtq is the initial time. Since 
| T;| 3> 1, for r > | Ti| the so-called weak- coupling asymp- 
totic approximation may be applied—. The asymptotic 
of the transition probability is 



P T ~ 1 - e" 



-e-™ 2 /Vl - e"™ 2 cos^(r), 
r 



(10) 



where 



cos^(r) = - +t 2 + — r 2 ln2 



, or 
argT( 




Adiabatic approximation 



-30 -10 10 20 30 40 
T 

FIG. 4: (Color online) Probability of transition |0) — > |1) as 
function of dimensionless time r (u — 3). Red line: proba- 
bility of adiabatic transition, P a dij). Blue line: probability 
of transition P T obtained from the exact solution of the LZ 
problem. 



In the adiabatic approximation the probability of the 
system to remain in the ground state may be described 
by 



Pad{r) = \(u-{t 



1 



Vr 2 + a; 2 



(12) 



The condition for adiabatic evolution, required by the 
adiabatic theorem, is u> 2 3> 1—. The probability of 
remaining in the ground state at the end of evolution 
(t — > oo) is given by P a d = |C'i(+oo)| 2 . For slow evo- 
lution we can use the LZ formul a 21 ' 22 to describe the 
probability of adiabatic evolution: 



Pad = 1 - e 



(13) 



In Figs. 01 [5] the probability of the adiabatic transition 
(red line) and the results of the exact solutions (blue 
line) are depcted. As can be seen for u> — 3 there is a 
good agreement between exact solution and the adiabatic 
formula (fT2"]l . For u> = 0.5 the disagreement takes place. 

Using the so-called adiabatic-impulse (AI) 
approximatio n 13 ' 15 , qualitatively, the dynamics of the 
Landau-Zener model can be described by the Kibblc- 
Zurek theory of nonequilibrium phase transitions^—. 
The Al-approximation assumes that the whole evolution 
can be divided in three parts and up to the phase factor 
the wave \u(t)) function approximately can be described 
as 



r € [— oo, — t\ 

t e [-f,f] 

r e [f , +ool 



|«(t)> « |«_(r)) 
|«(r)) « |u_(-f)) 
|(u(r)|u_(r))| 2 = const 



where the time f, introduced by Zurekl, is called the 
freeze-out time and define the instant when behaviour 
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FIG. 5: (Color online) Probability of transition |0) — \ |1) as 
function of r (tj = 0.5). Red (upper) line: probability of 
adiabatic transition. Blue (lower) line: transition probability 
obtained from exact solution of the LZ problem. 



of the system changes from the adiabatic regime to an 
impulse one where its state is effectively frozen and then 
back from the impulse regime to the adiabatic one. 

If the evolution starts at moment 7$ -C — f from 
the ground state, the equation for determining f reads 
ttt/2 — l/gap(f) (for details of calculation see ReL^ 3 -), 
and its solution is given by 




(14) 



Using the relation r = tq(P/P c — 1), we find that the 
change of the adiabatic regime to a non-adiabatic one 
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FIG. 6: (Color online) Probability, P e , of finding the system 
in the excited state (a; = ttuj 2 ): blue line - P e m Pai (I16|l . 

_ 2 

dashed red line - the LZ expression, P e — e nw . 



occurs when the pressure is Pi = P c (l — t/tq), and the 
non-adiabatic evolution becomes the adiabatic evolution 
again when the pressure increases up to P2 = P c (l + 
t/to). From here we find that within the interval of 
pressure, AP = P 2 — Pi = 2F c f/r , the behaviour of the 
system is described by the impulse regime. 

Employing Eq. (|14p . we approximate AP for fast 
(uj 2 <C 1) and slow (uj 2 3> 1) transitions as 



1 



AP 

~p7 



uj z < 1 



/7TT 
' w 2 » 1 



(15) 



TTUJTq 



In the AI approximation the probability, P e , of find- 
ing the system in the excited state at t/ S> f can be 
calculated as follow o 13 i 15 : 



P e ttP AI = \( U+ (T)\U-(-t))\ 2 = 

Substituting f from fl"4")) , we obtain 



Pai 



c 2 + x\Jx 2 + 4 + 2 



(16) 



(17) 



where x — uuj 2 . 

For uj 2 < 1, from Eq. (JTTJ) it follows P A i «1- vrw 2 . 

In the first order this coincides with the result predicted 

2 

by exact LZ formula: P e = e" 1 ™ . For the adiabatic 
evolution, uj 2 3> 1, we obtain Pai ~ l/7r 2 a; 4 (See Fig. 
[HI). As can be seen, the AI approximation is good enough 
for uj 2 < 1 and in the limit u> 2 ^> 1. 

Comparison with experimental data. - To com- 
pare our theoretical finding with the experiments on spin 
crossover under the high pressure we use the data from 
Refsi - — . The typical value of the critical pressure is 
P c = 50GPa, and the rest of parameters are taken as 
follows: £ = eV, p = 0.01 eV and tq w 10 4 s. The 
computation yields: tq sa 10 9 and to ~ 10 7 . Thus the 
spin crossover under slowly changed pressure realized 



in the cited experiments is a highly adiabatic process 
(cj 2 1). Using (fl~5"l) , we find the that the domain of 
non-adiabaticity is defined by AP w 10 _7 Pa. The cor- 
responding interval of time t « l/(7rpro) is t rj 10 -21 s. 



B. Quench dynamics under shock-wave load 

In this section we study the quench dynamics in the 
spin system under time-dependent pressure, P(t). We 
assume that at the initial moment of time P(U) — 0, and 
at the end of evolution P(t/) = Po- Further it is conve- 
nient to present P(t) as P(t) = Pos(t), where s(ti) = 
and s(tf) = 1. 

One can observe that the first term in the Hamiltonian 
([2]), yielding contribution to the total phase factor of the 
the wave function, does not affects the dynamics of the 
system and may be omitted. Setting for simplicity tp = 
and using Eq. (J3]), one can recast the time-dependent 
driving Hamiltonian as follows: 



Ur{t) = 



£ (1 - as(t)) 
P 



-e (l - as(t)) 



(18) 



where a = Pq/P c . For a given a, the crossover occurs at 
the critical point s c — s(t c ) defined as s c = I /a. 

In what follows we specify the pressure as a pulse with 
the shape determined by 



P = 



0, t < 

P tanh(crf), t > 0, 



(19) 



where Po is the pulse height. For Po = P c , expanding 
P(t) near of the critical point up to the first order, we 
obtain the related LZ problem: P(t) = P c (l + at). From 
here we obtain the corresponding LZ adiabaticity param- 
eter as ujlz = p/\/ct£o- 

Applying the adiabatic theorem we find that condition 
for adiabatic evolution can be written as follows: 



1 

— max 

p 1 0<P<P 



< 1. 



Next using Eq. (fl9|) . we obtain 



0<s<l 



oj9(l 



■p 2 



(20) 



(21) 



where /3 = p/so- In terms of the function w(s) defined as 



1 



a/?(l 



^ 2 v/(l-as) 2 +/3 2 ' 



(22) 



the condition of adiabaticity (|2"Tj) can be recast to the 
following: 



max 



1 



0<s<l \uj 2 {s) 



<. 1 



(23) 



6 





FIG. 7: (Color online) Parameter of adiabaticity as function 
of s. From the right to the left: a = 10, 5, 2.5, 1.25 (iv LZ = 1). 



t(ns) 

FIG. 9: (Color online) Dependence of transition probability 
P T on time t, for different initial conditions (Pq — P c ). From 
up to down: |Co(0)j 2 = 0, 0.25, 0.5, 0.75, 1 (a = 10ns" 1 ). 
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FIG. 8: (Color online) Dependence of /(a) on the amplitude 
a = P /Pc (0 = 0.01). 



FIG. 10: (Color online) Transition probability P T for different 
values of the pulse height as function of time t. From up to 
down: a = P /P c = 1.01, 1.001, 0.999, 0.99 (a = 10 ns" 1 ). 



In Fig. [7] the function 1/uj 2 (s) is depicted {ojlz =!)• 
We observe that 1/lu 2 has a maximum in the critical 
point. This is in agreement with the general observation 
on breaking-down of adiabaticity in the neighbourhood 
of the critical point. 

For j3 <C 1 we obtain the following estimate 



/(«) 



(24) 



where 



2a/3VT' 



(i + vi _ 7) v / T 7 fl 7 +F 

a — 1/a, a > 1 



a < 1 



(25) 



In Fig. [8] the function /(a) is depicted. We observe that 
/(a) -C 1 for a < 1 and it is increasing function, f(a) > 1 
for a > 1. 

In Figs. H] - [TT] we present the results of numerical 
calculations for different choice of parameter a. In all 
calculations the parameters were chosen as follows: Bq = 



leV, p = O.OleV and a = 10 ns" 1 
and for j3 — p/so, we obtain j3 



This yields lolz = 1, 
0.01. 



Fig. [S] shows that for Po = -Po independently of the 
initial conditions, the final state of the system is defined 
by the equal mixture of LH and LS states. The reason 
for this phenomena is in growing quantum fluctuations 
which lead to the mixture of LH and LS states at the 
critical pressure P c . In Figs. [l"0l and [TT1 the dependence 
of transition probability on time and pressure, respec- 
tively, is depicted for different values of the dimension- 
less amplitude of the pulse, a = Pq/P c . In the interval 
of pressure, < P < 1.25P C , the evolution is highly adi- 
abatic, and for Po ~ L25P C the system comes to the LS 
state with the minimal density of defects. However with 
the increasing of pressure amplitude, Pq, the probability 
to pass to the LS state is decreasing, and the system re- 
mains in the mixture of the HS ans LS states. Further 
increasing of the amplitude leads to the domination of 
the HS population over LS population. 
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FIG. 11: (Color online) Transition probability P T as func- 
tion of scaled pressure s for different values of the pulse 
height. From up to down: a = P /P c = 1.25,2.5,5,10,20 
(a = 10ns" 1 ). 



IV. SUMMARY AND CONCLUSIONS 

We have analytically studied the properties of spin 
crossover under the high pressure. We showed that at 
static loading and P = P c occupation numbers of both 
HS and LS states are equal at zero temperature (T = 0) , 
Jigs = n LS — 0.5. For P < P c we obtain Uhs = 1 an d 
tils = 0, while for P > P c one has nns = and tils = 1. 
Static transition at T = is a sharp quantum phase tran- 
sition with geometric Berry-like phase being the order 
parameter—. Finite temperature removes singularity in 



the tihs(P) dependence. Thermal fluctuations between 
two states |0) and |1) results in a smooth crossover in- 
stead of the quantum phase transition at T = 0. 

To verify our theoretical predictions we have performed 
numerical simulations. Fig. [10] shows that the tempo- 
ral quantum fluctuations have an effect similar to the 
thermal fluctuations. Small deviations of the shock wave 
amplitude Pq/P c from the unity results not in a sharp 
change of the probability P T either to zero or to unity 
but to a continuous deviation of the P T from the 0.5 
value. At the same time for Po/P c = 1 any initial dis- 
tribution of the HS and LS states will end its evolution 
in the equilibrium state njjs — n LS — 0.5 (Figj9]). This 
conclusion is valid for any choice of parameter a. 

Another dynamical effect we would like to discuss is the 
adiabaticity violations near the critical pressure (Fig lTTI) . 
This is a manifestation of the general Kibble-Zurek the- 
ory. The results shown in Fig llll are contra intuitive at 
the first glimpse. The shock wave with larger amplitude 
has smaller final probability for spin crossover and larger 
probability to stay in the initial HS state. To understand 
this effect one should note that the characteristic scale of 
the pressure increase is given by the factor aa for the 
wave with amplitude a. Thus larger amplitude wave also 
is faster. 
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